The critical behaviour of the Ising model on the 4-dimensional lattice 
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In this paper we investigate the nature of the singularity of the Ising model of the 4-dimensional 
cubic lattice. It is rigorously known that the specific heat has critical exponent a = but a non- 
rigorous field-theory argument predicts an unbounded specific heat with a logarithmic singularity 
at To. 

We find that within the given accuracy the canonical ensemble data is consistent both with a 
logarithmic singularity and a bounded specific heat, but that the micro-canonical ensemble lends 
stronger support to a bounded specific heat. 

Our conclusion is that either much larger system sizes are needed for Monte Carlo studies of this 
model in four dimensions or the field theory prediction of a logarithmic singularity is wrong. 
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I. INTRODUCTION 

In dimension D > 5 it is known from [1, 2] that the 
Ising model on the cubic lattice exhibits mean-field crit- 
ical exponents at the critical temperature. Even earlier 
it was shown [3] that the specific heat obeys the mean- 
field exponent a — for D > 4, and that for D > 5 the 
specific heat is in fact bounded at the critical point. For 
D = 4 the rigorous results which determine that a = 
are not strong enough to show that the specific heat is 
bounded. In fact methods from field and renormaliza- 
tion theory predict that the specific heat should diverge 
as (In \T — Tel) 1 / 3 but this has not been possible to prove 
rigorously. There are thus, at least, two possibilities here, 
either the specific heat is bounded in D — 4 as well or it 
diverges logarithmically. 

Earlier studies of the critical behaviour in 4-dimensions 
include [4-6], using Monte Carlo methods, and [7], using 
series expansion and extrapolation. There has also been 
some recent controversy [8-10] regarding the consistency 
of field theoretical predictions and Monte Carlo data. 

Using a standard Monte Carlo approach to detect a 
divergence of the form (In IT-Td) 1 / 3 is difficult since the 
quantity will remain quite small for a large range of the 
lattice size L, thereby making it difficult to use sampled 
data to clearly distinguish between different asymptotic 
behaviours. 

In an attempt to get around this problem we have in- 
stead studied the microcanonical density of states of the 
model, following the methods used in e.g. [11, 12]. The 
finite-size effects of the canonical ensemble have two com- 
ponents; that coming from the fact that only a certain 
discrete set of energies are available in finite discrete sys- 
tems, and that coming from finite-size effects of the den- 
sity of states. The microcanonical ensemble is affected 



* Email: phl9kth.se 

t Email: klas . markstrom9math . umu . se 



by only the latter effect. 

A divergence in the specific heat means that the second 
derivative of the density of states must become at the 
critical point. The surprising simulation result is that 
this value is in fact increasing with the lattice size at 
the critical point and the best fit to the data is that 
it converges to a non-zero value, thereby also giving a 
bounded specific heat in the limit. 

In order to make sure that this was not an artifact 
caused by our simulation software we wrote two separate 
programs, one for the Metropolis algorithm and one using 
the Wolff-cluster algorithm [13], to sample at interleav- 
ing lattice sizes, but no systematic differences could be 
seen. We also tried to push the simulations to large lat- 
tices, reaching L — 80. Our simulations give estimates for 
the critical exponents which agree well with the rigorous 
mean-field values and a value for the critical temperature 
which agrees well with earlier studies. 

Hence our conclusion is that either lattice sizes larger 
than L = 80 are needed to see the asymptotic behaviour 
of the specific heat or the specific heat is in fact bounded 
at the critical point. Finding a way to settle this issue is 
of prime importance since, as discussed in e.g. [10], this 
would have consequences for the renormalization tech- 
niques used to bound the Higgs mass. 



II. NOTATION AND BASIC DEFINITIONS 

The lattice studied here is the cartesian graph prod- 
uct of four L-cycles, that is, an L x L x L x L-lattice 
with periodic boundary conditions on n = i 4 vertices 
and m = 4I/ 4 edges. We have collected sampled data 
using the sampling scheme described in [11] for linear or- 
ders: L = 4,6, 8, 10, 12, 16, 20, 24, 32, 40, 48, 56, 60, 64, 80. 
For most orders we used the Metropolis single-spin flip 
method with measurements of local energies after every 
sweep. Since the flip-rate near the critical temperature 
is about 63% there will be no strong dependency be- 
tween measurements of local energies. For comparison 
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we also employed the Wolff-cluster method for the cases 
L = 10,20,40,60, flipping clusters until an expected L 4 
spins were flipped. 

The energy E of a state a = (01, . ..,<j n ), with er^ = 
±1, is defined as E(a) = ^ <t% <Jj , with the sum taken 
over all the edges {i,j}, and the magnetisation M is de- 
fined as M(a) = <Ji with the sum taken over all the 
vertices. 

We have two classes of quantities. First the combina- 
torial quantities from the microcanonical ensemble which 
depend on the energy U . Especially the coupling K 
is of interest here, defined as K(U) = —dS/dU where 
S(U) = {\na{E))/n for U = E/n and a{E) denotes the 
number of states a at energy E. How to obtain the cou- 
pling from sampled data is described in detail in [11] and 
error estimation in [14]. 

The canonical, or physical, quantities are obtained as 
cumulants, or derivatives of In Z(K, H) with respect to 
K or H (the external field), where Z is the partition 
function. All quantities are measured with the external 
field switched off, ie H — after the relevant derivative 
is taken. 

At this point we introduce the notation Cj — 
((A — (A)) 1 ) for the ith central moment of a random 
variable A, where (A) is the mean value. The kth cu- 
mulant of E is then the fcth derivative of In Z with respect 
to K, where Z is the partition function. Recall that the 
first cumulant is (A), the second is 02(A) = Var(A), 
the third is 03(A) and the fourth is 04(A) — Sc^A). 
The internal energy is then U{K) = (E) jn and the spe- 
cific heat is C(K) = Var (E) jn. Note also that the 
susceptibility \ = Var(M)/n = (A/ 2 ) jn has no lo- 
cal maximum, whereas the (spontaneous) susceptibility 
X = Var (|M|) jn does. Analogously we define the mag- 
netisation as [i = (M) jn and the spontaneous magneti- 
sation as Ji — (\M\) jn. 

III. PHYSICAL QUANTITIES 

Let us begin by showing some plots of a few physical 
quantities near the critical coupling. Figure 1 shows the 
magnetisation ji{K). In Figure 2 we show the specific 
heat C(K) for several lattice sizes. 

A. Critical points and exponents 

First we establish a high-precision estimate of the crit- 
ical coupling K c . This is done by determining the critical 
points for a number of different quantities, listed below, 
for each system size. The critical points in question are, 
with one exception, the locations of various maxima or 
minima of eg cumulants. To these points we fitted a 
simple scaling law of the form Co + ci L~ A . By selecting 
points for L > L m ; n for different L m ; n we can then obtain 
several (for L m i n = 24, 32, 40, with a few exceptions) dif- 
ferent estimates of the fitting parameters. As a rule we 
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FIG. 1: (Colour online) Magnetisation /j.(K) for lattice sizes 
L = 6, 8, 10, 12, 16, 20, 24, 32, 40, 48, 56, 60, 64, 80. 
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FIG. 2: (Colour online) Specific heat C(K) for lattice sizes 
L = 6, 8, 10, 12, 16, 20, 24, 32, 40, 48, 56, 60, 64, 80. 



received very good fits deeming a higher order correction 
term unnecessary. The sought parameter is of course Co- 
Taking the median of these gives a final estimate of K c 
for that particular quantity. Repeating this for all quan- 
tities, a grand total of 15, allows us to make a statistical 
analysis of them. We have used the median as the esti- 
mate, with the first and third quartile as error estimates. 
In short, we take the median of the medians, very much 
like in [12]. 

The points scale very nicely with the linear order using 
only the simple expression above, see Figure 3. The re- 
sulting estimate is K c = 0.1496947±5-10~ 7 . This falls in- 
side the by now rather old estimate K c = 0.14965±5-10 -5 
found in [16] and agrees with the estimate from [4]. 

The critical points in question are the locations of 
the following; the maximum of the specific heat C and 
susceptibility x, maximum and minimum of the cu- 
mulants c 3 (E)/n, c 3 (\M\)/n, (c 4 (E) - 3c 2 (E))/n and 
(c 4 (|M|)-3c|(|M|))/n, maximum of dft/dK, d\njl/dK, 
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d\ax/dK and dQ/dK, where Q is the Binder cumulant 
l-(Af 4 )/3 (M 2 ) and finally the crossing point between 
Ql and Ql/2- See eg [15] for a discussion of the last four 
quantities. 
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FIG. 3: (Colour online) The critical points vs 1/L with fitted 
curves. 

The expression above also provides us with estimates 
of the exponent v. The location of a critical point K* L 
should deviate from K c as roughly K* L — K c ~ L~ 1//iy , 
again see [15]. Repeating the median-of-the- medians ap- 
proach gives A = 1/v = 2.00 ±0.03 where the bounds are 
again based on the 1st and 3rd quartile, thus rendering 
us v = 0.50 ± 0.01. The Josephson inequality tells us 
that a > 2 — Dv, and hence our midpoint estimate gives 
a > ± 0.04 for D = 4, since a = [3] our data is in 
good agreement with the rigorous results. Similarly an 
estimate of j3 — 0.50 is found, and the mean- field value 
is/3=|. 

Having established an estimate of K c we can now esti- 
mate the internal energy Ul(K c ) and again fit the scaling 
formula above to these data for L m i n = 24,32,40. The 
different Co, and thus the asymptotic values of U c , end 
up inside the interval 0.77053 ± 4-10~ 5 . 



B. Critical values 

Our aim is now to try to distinguish between the two 
possible scenarios, either we have a logarithmic singular- 
ity or the specific heat is bounded at T c . We attempt to 
do this by making least-squares fits to the data for two 
different forms of the fitting function. 

According to scaling theory, see [17], the maximum 
specific heat C max is proportional to v'ln L. For L > 
12 this seems plausible given our data. In Figure 4 we 
show C max versus v^ln L together with a fitted straight 
line, y = 115 x — 56.7, and indeed they line up rather 
convincingly. The reader should note that C max grows 
very slowly indeed. 

For the bounded scenario we try a fit where C max is 
proportional to a power of L. A least-squares fit of both 



constant and exponent gives 150.49+ 180. 5L -0 496 . We 
show this in the inset of Figure 4. The fact that the ex- 
ponent is negative would of course mean that the specific 
heat is finite in the limit. 

For both models there is some variation in the coeffi- 
cients and the exponent if one makes the fit to different 
subsets of the data points, but no drastic changes. An 
attempt with evaluating the specific heat and the suscep- 
tibility at the asymptotic K c for each linear size instead 
gave a very similar behaviour to that of their maximum 
value. 

To the eye both fitting functions work reasonably well 
and we simply find that the canonical ensemble data can 
not strongly distinguish the two scenarios. 
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FIG. 4: (Colour online) C max 1 
together with the fitted curves. 
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IV. COMBINATORIAL QUANTITIES 

With regards to the microcanonical ensemble the two 
scenarios will be that either K'{U) goes to at U c or it 
converges to a finite positive value 

Figure 5 shows the microcanonical quantity K(U) and 
in Figure 6 its derivative is shown, both together with 
zoomed- in versions near the critical energy U c . Most of 
the sampling was done for energies close to the critical 
one for the given value of L so the curves become noisier 
further away from U c . 

The minima do not at all seem to approach zero as 
they do for d = 2 [18] and d = 3 [12]. In fact the be- 
haviour here is qualitatively different in that the values 
are actually increasing rather than decreasing. 

It is known, see e.g. [11], that the specific heat corre- 
sponds to 1/K'(U). Thus limu^. Uc K'(U) = if and only 
if ]im.K^K a C{K) = oo. Figure 7 shows the minima ver- 
sus 1/L together with a fitted line y = 0.00419-0.0151 x, 
suggesting that the minimum approaches a maximum 
0.00419. 



The optimal exponent of 1/L, naturally, depends to 
some extent on which data points are used. Using a least- 
squares fit to different subsets of the data for L > 6 
gives exponents between (roughly) 0.9 and f.5. More 
specifically, if we check all subsets of the data with L > 6 
on between 10 and 12 points a median exponent of 1.25 
is received and for cq the median value was 0.00406 with 
first and third quartiles 0.00402 and 0.00413 respectively. 
The extremal values for Co are 0.0038 and 0.0044. If we 
instead use all the data points for L > 8 we obtain the 
exponent A = 1.147 and c = 0.00411. 



for the given lattice sizes the canonical ensemble can not 
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FIG. 5: (Colour online) Coupling K(U) for L > 6. 




FIG. 6: (Colour online) Coupling K'(U) for L > 6. 



FIG. 7: (Colour online) The minimum of K'(U) vs 1/L, to- 
gether with the fitted curve. 



conclusively distinguish between the two scenarios, and 
in an attempt to circumvent this we have instead turned 
to the microcanonical ensemble. 

There are two reasons for why the microcanonical en- 
semble could give clearer results in this situation, the first 
predicted and the second unexpected. 

First, the canonical ensemble is expected to have larger 
finite size effects than the microcanonical ensemble. To 
see this we may consider an idealised example where, for 
a finite system, S(U) at each energy U is identical to 
the limit as n — > oo. Here the density of states has no 
finite size effects at all, apart from only being defined at 
certain discrete set of values of U . However because of 
the discrete energies there will still be finite-size effects 
in the corresponding canonical ensemble. 

Secondly, a divergent specific heat means that K'(U) 
goes to at U c , and as we have found the minimum value 
of K'(U) is actually increasing rather than decreasing. 
This gives us a qualitative signal, rather than a weak 
quantitative one, that the specific heat actually converges 
to a finite value. 

Our conclusion is that either much larger systems are 
needed to see the asymptotic behaviour of this model, 
and this possibility can only be ruled out by a rigorous 
convergence result, or the specific heat is in fact bounded 
at U c , thus contradicting the renormalization theory pre- 
diction. 



V. CONCLUSIONS 

We have studied the two proposed scenarios for the 
critical behaviour of the specific heat of the 4-dimensional 
Ising model. This has been done in both the canonical 
and the microcanonical ensembles. We have found that 



Acknowledgements 

This research was conducted using the resources of 
High Performance Computing Center North (HPC2N) 
and the Center for Parallel Computers (PDC). Thanks 
are also due to the referees for their constructive criti- 
cisms. 



■5 



[1] M. Aizenman, Phys. Rev. Lett. 47, 1 (1981), ISSN 0031- 
9007. 

[2] M. Aizenman, Comm. Math. Phys. 86, 1 (1982), ISSN 
0010-3616. 

[3] A. D. Sokal, Phys. Lett. A 71, 451 (1979). 

[4] E. Bittner, W. Janke, and H. Markum, Phys. Rev. D (3) 

66, 024008, 8 (2002), ISSN 0556-2821. 
[5] H. W. J. Blote and R. H. Swendsen, Phys. Rev. B 22, 

4481 (1980). 

[6] E. Sanchez- Velasco, Journal of Physics A: Mathematical 
and General 20, 5033 (1987), URL http://stacks.iop. 
org/0305-4470/20/5033. 

[7] M. Hellmund and W. Janke, Physical Review B (Con- 
densed Matter and Materials Physics) 74, 144201 
(pages 9) (2006), URL http://link.aps. org/ abstract/ 
PRB/v74/e 144201. 

[8] P. Cea, M. Consoli, and L. Cosmai (2005), hep- 
lat/0501013. 

[9] P. M. Stevenson, Nuclear Phys. B 729, 542 (2005), ISSN 
0550-3213. 

[10] J. Balog, F. Niedermayer, and P. Weisz, Nuclear Phys. 



B 741, 390 (2006), ISSN 0550-3213. 
[11] R. Haggkvist, A. Rosengren, D. Andren, P. Kundrotas, 
P. H. Lundow, and K. Markstrom, J. Statist. Phys. 114, 
455 (2004). 

[12] R. Haggkvist, A. Rosengren, P. H. Lundow, K. Mark- 
strom, D. Andren, and P. Kundrotas, Adv. Phys. 56, 
653 (2007). 

[13] U. Wolff, Phys. Rev. Lett 62, 361 (1989). 
[14] P. H. Lundow and K. Markstrom, arXivxond- 
mat/0612465. 

[15] A. M. Ferrenberg and D. P. Landau, Phys. Rev. B 44, 
5081 (1991). 

[16] D. S. Gaunt, M. F. Sykes, and S. McKenzie, J. Phys. A 
12, 871 (1979). 

[17] E. Brezin, J. C. Le Guillou, and J. Zinn- Justin, in Phase 

transitions and critical phenomena, Vol. 6 (Academic 

Press, London, 1976), pp. 125-247. 
[18] R. Haggkvist, A. Rosengren, D. Andren, P. Kundrotas, 

P. H. Lundow, and K. Markstrom, Phys. Rev. E 69, 

046104 (2004). 



